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^) • The Standard Model constraints on a which can be derived from the B — > vrvr decays are revis- 

fH , ited in some depth. As experimental inputs, the B^ — > vr+vr", B^ — > vr+vr'^ decays complemented 

Qh! by the B^ — > vr'^vr'^ decays, the CP parameters S'tttt and Ctttt, and/or the value of a as determined 

Q^' by the global CKM fit are used. The constraints discussed here are model independent in the 

sense that they rely only on Isospin symmetry, following the Gronau-London proposal. A new 
bound on B^^ and the function Coo{B^^) are introduced. While another bound applied to BABAR 
results is shown to imply that cos(2aeff) is negative. The Grossman- Quinn bound is rediscussed. 
A close form expression is given for a as a function of the measurements. Various scenarii for 
C^ ' the future of the isospin analysis are explored. To probe the Standard Model the {B^^, Cqo) plane 

is introduced. 
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1 Introduction 



The purpose of this note is to revisit the Standard Model constraints which can be derived 
from (and on) the B ^ tttt decays following the Gronau-London proposal [1]. Much have 
been written already on this topic, for example [2, 3, 4], but fresh developments are 
presented here. 

This study leads to correct a strong bias: the B ^ titt isospin analysis is not doomed 
neither to wait for ultra high luminosity, nor to have a small branching ratio B^^ (corre- 
sponding to the B^ -^ TT^TT^ decay) to be able to constrain a. The current period of time 
is very exciting: anything can happen in the next three years, including a rather precise 
determination of a. 

In Section 2 is recalled the present experimental status: the values of the three branch- 
ing ratios {B+- = B{B^ -^ tt+tt-), B+^ = B{B+ -^ n+n^) and i3°°) and the two time- 
dependent asymmetry parameters (S'^tt and Ctttt)- The statistical meaning of the BABAR 
and Belle results for S'tttt and (7^,^ is discussed, considering two different statistical treat- 
ments. 

In Section 3, the plague of the mirror solutions is explicited. Section 4 then turns to 
some algebra. A new bound on B^^ is presented Eq.(24), from which the sign of cos(2Q;cff) 
is inferred, a function Coo{B^^) is explicited Eq.(35), a close form expression for a as a 
function of the measurements is given, and the meaning of the Grossman-Quinn bound 
on I a — tteflf I is discussed. 

In Section 5 the perspectives for measuring a are pondered. First an exploration of 
the space of still open possibilities is attempted. Then two scenarii for the high luminosity 
regime (10 ab~ ) are considered. For the sake of completeness a nasty and a nice scenario 
are considered, with 500 fb^^. 

The analysis is reconsidered from another point of view in Section 6. Instead of 
the usual p and fj variables, the branching ratio i3°° and the asymmetry term Cqq are 
taken as key variables to probe the Standard Model, when dealing with B ^ ttti decays. 
The present experimental status and the high luminosity potential are discussed in this 
context. 

The Electro Weak Penguin complications are neglected throughout the note, having in 
mind that, when accounted for, they should not utterly blur the picture. 

The effects of radiative corrections [5, 6], which imply a 5 to 10% increase on the 
published branching ratios are not taken into account in the present analysis. 



2 Comparison of BABAR and Belle results 

The two detectors BABAR and Belle located respectively at the asymmetric B Factory 
PEP-II and KEK-B measure branching ratios and time- dependent asymmetries for the 
B — > TTTT decays. If the BABAR values, for the time-dependent asymmetries, are in the 
physical region. Belle ones are not. A special statistical treatment is then advised to deal 
with such a situation where a measurement lies beyond the physical boundaries: this is the 
purpose of this section. First are recalled in Section 2.1 the values measured by BABAR 
and Belle. The world average (W.A.) ones are given as well [7]. Section 2.2 discusses the 
results using the standard statistical treatment. Section 2.3 turns to the boundary effect 
via a particular treatment. 

2.1 Present experimental status 

The most recent values of BABAR and Belle are summarized in Table 1. All branching 
ratios are quoted in the 10~^ unit, as in the whole document. For the world average values 
of the two time-dependent asymmetries terms S^tt and C^tt, the discrepancy between 
BABAR and Belle is not taken into account and the results are merged together at face 
values. That is to say the raw measurements are taken, without accounting for the fact 
that physical values should be in the disk S^^ + C^^ < 1. Since Belle results are outside 
the physical disk, in some instances modified and fictitious results are used instead, where 
the Belle measurements are forced back onto the unity circle: they are denoted Belle* in 
the present document. 





B+- (10-6) 


i3+o (10-'^) 


^00 (10-6) 


L^TTTT 


^TTTT 


P 


BABAR 


4.7 ±0.6 


5.5 ±1.2 


2.1 ±0.7 


-0.19 ±0.20 


-0.40 ±0.22 


-0.02 


Belle 


4.4 ±0.7 


5.0 ±1.3 


1.7 ±0.6 


-0.58 ±0.17 


-1.00 ±0.22 


-0.29 


Belle* 


- 


- 


- 


-0.55 ±0.17 


-0.84 ±0.22 


- 


W.A. 


4.55 ± 0.44 


5.18 ±0.77 


1.9 ±0.5 


-0.46 ±0.13 


-0.73 ±0.16 


-0.17 



Table 1: Input values to the fits for the branching ratios [8, 9, 10, 11] and for the CP 
violating asymmetries [12, 13]. The notation Belle* refers to the Belle results for the 
branching ratios, but where the measured values for St^tt and Ct^t^ are replaced by the 
values quoted by Belle for S^^^^ and C^'^^^ (Section 2.3). The world average branching 
ratios are obtained using in addition CLEO results [I4]. 



2.2 The unbounded (standard) statistical treatment 

A proper method to compute a confidence level (CL) associated to a given theoretical 
hypothesis (hyp) and an experimental result (exp) is the following: one defines a test 
statistics w(hyp,exp) which bears the property that the larger its value, the smaller the 
likelihood of the hypothesis (hyp) to be true, given the experimental result (exp). One 



denotes /hyp(w) the probability density function of the variable w, the hypothesis being 
assumed to be true. The confidence level is then computed as: 



CL(hyp,exp) = / fi,yp{w)dw . (1) 

Jw(hyp,cxp) 

When the theoretical hypothesis is unbounded, a standard treatment is to identify w with 
the x^ measuring the agreement between theory and experiment. When the theoretical 
hypothesis is bounded, this standard treatment is still valid. However, it suffers from 
a drawback: if, because of a statistical fluctuation, the experimental result lies far away 
from the allowed region, the confldence level might be very small, whatever the hypothesis. 
The underlying theory is therefore ruled out. The drawback stems from the fact that one 
may wish to set a confidence level assuming that the underlying theory is true. That is 
to say, one does not question the vahdity of the theory as a whole, but one is interested 
to reach a statement on some of its theoretical parameters. 

2.3 The bounded statistical treatment 

The standard, x^-based, statistical analysis of the S',r7r and Ctttt measurements may lead 
to the exclusion of the whole physical domain S^^ + C^^ < 1 if statistical fluctuations 
bring the measurements far enough outside it. Because Belle values enter deep into the 
unphysical region, the hypothesis of no CP violation in the B ^ mi {Stttt = Ctttt = 0) is 
ruled out at 2.48 x 10^^% CL by her analysis [13]. Here is presented a statistical treatment 
which by construction cannot rule out the whole physical domain. It incorporates the 
(mild) assumption that at least Quantum Mechanics is right and hence that S^^+C^^ < 1. 
Let us denote: 

• S'^^P and C^^P , the theoretical values one would like to set a CL on. They are within 
the physical domain in the {S.„n, Ctttt) plane, 

• S^"^ and C°^P , the measured values. They can be anywhere in the {S^tt, Ctttt) plane, 

• S^^^^ and C^^^* , the theoretical values which are in the best agreement with the above 
measurements. They are within the physical domain in the (S'tttt, Ctttt) plane. These 
values are the ones used to replace S!^^ and C^"^, when considering the Belle* results, 

• x^(Gxp|hyp) is the x^ o^is gets when comparing the measured values with the hy- 
potheses S^IP and C^y/, 

• x^(exp|best) is the minimal x^(exp|hyp) one gets when exploring the physical do- 
main. The corresponding values of ^tttt and Ctttt are the above deflned S^'^^ and 

/^best 

^TTTT ) 

• w = x^(exp|hyp) — x^(exp|best) is the test statistics used to deflned the Confldence 
Level in the presently advocated treatment. 

If S'^^ and C^^P are inside the physical domain then x^(exp|best) = 0. Whatever S^^ 
and C^^P are, by construction, there is a point in the physical domain where w = 0: this 
point is {S^l^^ ,C^^^^} . It is on the boundary of the physical domain if the measured values 
are outside ^. 



^The precise location of this point is not at the intersection of the unity-radius circle and the line 
joigning the measured point and the origin, but not far. 
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Figure 1: Confidence level at 32% and 1% corresponding to BABAR and Belle results, and 
their average. The calculation takes into account the fact that the theoretical values are in 
the physical domain S^^ + C^^ < 1. The ellipse located inside the shaded areas indicate 
the 32% bound obtained with the usual treatment. 

For each point {S'^^p,C^^p}, one computes w from the values of the measurement 
at hand. Then, one proceeds to a (large) toy Monte-Carlo simulation to compute the 
probability to get a w value larger than the one above. In principle, one should first 
obtain the pdf P(^^^p, Q^p | 8^1^,0]^^), which shape slightly depends on {S^yP^C^yP}, 
and Belle does that. In practice, a mere Gaussian approximation is used, which is fair 
enough for the present discussion. Applying this treatment which accounts for boundary 
effects, one gets the contours of Fig. 1 where are drawn the CL computed from BABAR 
results (blue). Belle results (red) and from their average. For each of the three above sets 
of values, three contours are shown: the 32% (shaded area) and 1% CL ellipses provided 
by the treatment, and the 32% CL (smaller) ellipse obtained with the standard treatment. 

BABAR results are still compatible with no-CP violation, but viewed from Belle or 
W.A. ones, the hypothesis of no-CP violation is excluded at a very high level. 



3 Isospin Constraints on B -^ rnv 

The constraints implied by SU(2) on the B ^ Tin decays [1] can be expressed as: 
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tt^tt'^ decays, respectively, 

• H, a and a are three unknown real (and positive) parameters which drive the strenght 
of the branching ratios, 

• ao, acs and 6 are three phases (ao is the true value of the angle a), 

• the phase convention chosen here is such that A'^^ is real positive. 

Although deff bears information on CP violation, its value is technically independent of 
the tto value: the presence of Penguins diagrams breaks down the relationship, if they 
are theoretically unknown, which is assumed here. In terms of the above six unknown 
parameters, the six physical observables take the form ^: 

-B+° = fx' (9) 

^00 ^ l(|^uu |2^|^" |2)^^2- ,2 + ^(a^ + a^)- V2(ac + ac)) (10) 
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where it proves convenient to use the two quantities c and c defined as: 

c = cos(0), = ao — 6 (14) 

c = cos(0), = cto + 5 — 2acflf • (15) 

The quantity 6*00, the analogous for 7r°7r° of S'^jr, could also be considered. Its expression 



is 



2 sm{2ao) + aa sm{2acfi) — o,\^sm{ao + S) — a^/2sin{ao — 5 + 2aefi) /-,„n 

°° " 2 + i(a2 + a2)-v^(ac + ac) ' ^ ^ 



^To take into account the fact that the hfetime of charged and neutral B are different, every branching 
ratio should be divided by the corresponding lifetime: B^ , ,B+° and B^^ become then respectively t^it: 

^^Tf and ^gu ■ AH numerical results are obtained with this modification, but for the sake of clarity the 
lifetimes do not appear explicitly in the formulae. The values used for the lifetimes are taken from [15]. 



However, its measurement (through Dahtz decays of at least one of the 7r°) imphes huge 
statistics and it is assumed that 6*00 is not measurable. The six remaining measurements 
B^^ , i3+°, i3*^°, S'tttt, C^tt and Coo can in principle fix the six unknowns, including a, at 
least up to multiple solutions. The phases a^s, S and ao enter the above six Eqs.(8-13) 
only through c, c and sin(2Q;eff), which are invariant through the three changes <-> —(f), 
(f) ^^ —0 and 2acs <-* tt — 2aes- Hence, there are 8 mirror solutions in the range [0, n] (16 
in the full range [0,27r]). These 8 mirror solutions are explicited in Table 2 (up to a ±7r 
shift to bring them back in the range [0, vr]) where the first row indicates the true solution 
for the triplet {a,S,aef[). 
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ftcflf 
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ftcflf 








2 


6 


ao 


«cflf 


-0 
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-ao + 2aeflf 


-5 + 2acff 


ftcflf 


-0 


-0 


4 


-6 + 2acff 


-ao + 2acflf 


ftcflf 





-0 


5 


f-«o 


f-^ 


f — ttcflf 


-0 


-0 
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f-^ 


f-«o 


f - «eflf 





-0 


^7 

1 


f + ao - 2acff 


f + 5 - 2aeff 


f -tteflf 








8 


f + (5 - 2aeff 


f + ao - 2acflf 


f — ttcflf 


-0 






Table 2: Mirror solutions for a in the range [0,7r]. All the 8 triplets (a,5,acs) yield the 
same values for c, c and sin(2aeflf). Solutions (5) to (8) are just 7r/2 minus solutions (1) 
to (4)- The 8 mirror solutions are strictly equivallent. They cannot be distinguished if no 
input is added to the ones of Eqs.(8-13), like Soo for example. 



It is worth remarking that the mirror solutions appear in pairs symmetrical around ac^: 



a(l) - aeff 
a (2) - aeflf 



aeflf - a(3) 
aeflf - a (4) 



(17) 
(18) 



and that, if the mirror solutions closest to the true value ao are solutions (2)-(3)-(4), 
denoting them a{ii,iii,iv) according to their closeness rank: 



a(n) — ao |<| a{iii) — ao |<| a{iv) — ao 



(19) 



It follows, from Table 2, that the fourth-rank mirror solution cannot be too close, nor too 
far away from ao = a{i), since it satisfies the double inequality: 



2 I a{ii) — ao |<| a{iv) — ao |< 2 | a{iii) — ao 



(20) 



Therefore, if the second-rank mirror solution a{ii) is near ao, while the third-rank mirror 
solution a{iii) is also creeping close, then the fourth-rank solution a{iv) is nor far either: 
in that case, ao will not be precisely measured. This topic is discussed in Section 4.4 in 
the light of the Grossman-Quinn bound. 



4 Algebra 

In this section one deals with some stark mathematics. The algebra involved being quite 
simple, only the results are presented: they are important for the next sections. It is 
demonstrated in passing that the sign of cos(2Q;eff) is firmly established by BABAR data, 
and the current CKM fit results. Below, are discussed the two types of bounds one may 
find in the literature: lower (and upper) bounds on i3°°, and upper bounds on | a — a^s \- 
These bounds are of some interest, but they can also be misleading. Furthermore, because 
of the limited precision on the measurements they use, all bounds (and the new bound 
of Eq.(24) is no exception) do not provide more than qualitative insights. To turn a 
bound value into a quantitative statement, one should append the CL it refers to, and 
this implies to have recourse to fits, which is done in Section 5. 

4.1 New bounds on B^^ 

New bounds on B^^ are presented, using the current knowledge on a. The value deter- 
mined by the standard CKM fit is 96 ±13° (1.67±0.23 rad) [16]. Using a presents several 
interests: 

• to help cornering i3°°. For the time being, the question is not yet to pinpoint a 
thanks to B ^ ttti decays, but to evaluate the chances one has to do so in the close 
or distant future. With that respect, one is longing for a robust estimate of i3°°, and 
every tools at hand should be used to secure it. The information on a must be used 
to this end as well, 

• to derive the sign of cos(2acflf) (Section 4.1.2), 

• to understand the behavior of the mirror solutions (Section 5.1). 

Assuming that a is known, the set of equations for the B ^ mr amplitudes implies the 
following lower bounds for the 5° -^ tt^tt^ branching ratio ^: 

?oo _ R+O 1 



B^^ = i3+^ + -i3+--^i3+oi3+-(l + 0±) (21) 



0± = Jl- Cl^ cos(2(ao - «cff)) (22) 



= ±V(l-sin^(2ao))(l-C2^-5^J+sin(2ao)5.. (23) 

where ± in 0± is the sign of the product cos(2Q;eflf) cos(2q;o). 

They come with the upper bounds B^^.^^ obtained by reversing the sign in front of the 
square root. However, the upper bounds are not very relevant, being much weaker than 
the direct experimental upper limit. 

4.1.1 The reBound 

The bound B^^.^j^ obtained when cos(2Q;efif) cos(2q;o) is taken to be positive is referred to 
as the reBound: 



They are obtained by deriving B with respect to 6 



0+- 




(24) 
(25) 


B^V = ^+° + ^B+- - ^B+^B+-{1 + 0+) 


= v/(l - sin2(2ao))(l - Cl^ - ^^J + sin(2ao)5.. 



The alternative bound i3^°_ is considered in the next section. 

The reBound can be compared with previous bounds [17, 18], the strongest lower one 
being the Gronau-London-Sinha-Sinha (GLSS) bound [19]: 



1 



^GLSS:i = ^^° + 2^^" - \/^+°^+- (1 + V 1 - Cl^) (26) 

which minimal value with respect to Ctttt is attained for Ctttt = 0: 

^S°LSS:i[min] = i3+°^-^/^j . (27) 

It comes in pair with the upper bound -Bo'Lgg.y obtained by reverting the sign in front of 
the square root. However, as obvious from Eq.(22) the reBound is more stringent than the 
GLSS bound, which it reproduces if no knowledge on a is assumed. The largest difference 
is obtained with the Belle* measurements (C^^ + S'^^ = 1), which, taken together with 
the value ao = 90°, imply the new lower limit B^^ > 2.3, while the weaker bound of 
Eq.(26) is i3°° > 0.7 (see Table 3). 

Using unadorned standard error propagation to estimate the precision of the bound, 
one gets for the reBound, a = 96°, and, 
BABAR values: 

s°V = o.6+°:l[^+°]l°:?[^+ll°:J[5..]t°:^[a.]tS:^[a] (28) 

Belle* values: 

s°V = i-9i:^[^+°]l^:J[^+-]l^:}[5..]lH[a.];HN (29) 

and world average values: 

^°V = i.iiJ[^+°]l^:}[^+1l^:^[^..]ii[a.]+^:^[«] (30) 

where the source of the uncertainty on the bound is specified within brackets. Whereas this 
tells that the main negative inaccuracy on the reBound value stems from the branching 
ratio B^^, one cannot infer much else: the errors are large and correlated, hence they 
cannot be added in quadrature. Anyway, one should keep in mind that the actual 90% 
CL lower limit on B^^ can hardly be guessed from back of the envelope calculations. 

The numerical results for the reBound for BABAR, Belle and world average are sum- 
marized in Table 3 for three different Oq values. In this table are also given numerical 
values of other bounds presented in the following sections. 

10 



Bounds 


"GLSSiI 


BZ+- rcBound 


i3°V: tough Bound 


^00/^+0 


'S'gq 


a 




90° 


96° 


110° 


90° 


96° 


110° 






BABAR 


0.6 


0.7 


0.6 


0.6 


5.9 


5.2 


3.8 


0.4 


0.4 


Belle* 


0.7 


2.3 


1.9 


1.2 


2.3 


1.9 


1.2 


0.4 


0.3 


W.A. 


0.7 


1.3 


1.1 


0.7 


3.8 


3.3 


2.2 


0.4 


0.4 



Table 3: Numerical values for the different bounds: the GLSS bound, the reBound, the 
tough Bound, the original (approximate) Grossman-Quinn bound and the exact Grossman- 
Quinn bound. Branching ratios are quoted in the 10^® unit. The last two values are 
discussed in Section 4-4- 

4.1.2 The tough Bound 

The hypothesis that the sign of cos(2Q;eflf) is opposite to the sign of cos(2q;o) is now made. 
A more stringent bound is then obtained: 






i3+° + ^B' 



B+OB+-(1 + 0^ 



(31) 



which is referred to as the tough Bound. For Belle*, since C^^ + S^^ = 1, the tough Bound 
is identical to the reBound. However, for BABAR and for the world average, reverting the 
sign makes a huge difference. The precision on the tough Bound for a = 96° is, for the 
BABAR values: 



B[ 



00 



for Belle*: 



B[ 



00 



;1- 






and for the world average: 



BZ- = S.Str.lB 



+0.5ri2+0l+0.0ri2+-l+0.4 



-0.0 



[Bn^ToiS: 



hO.2 



[a 



1+0.7 



TTTj-O.Sl'-^'r'rJ-O.? 



a\ 



(32) 
(33) 

(34) 



These tough Bound values are, for BABAR and W.A., firm enough to be barely consistent 
with the measurement of i3°° (Table 1): the hypothesis cos(2aeff) cos(2q;o) < is disfa- 
vored. On the other hand, the sign of cos(2q;) is well defined to be negative by the current 
state of the art of CKM fits. Hence, the tough Bound suggests cos(2Q;eff ) < 0. 

4.2 The Coo function 

The asymmetry term Coo (Eq.(13)) has not been yet measured due to the insufficient 
luminosity. From the knowledge of a, once the three branching ratios i3^~, iS^", B^^ and 
the two asymmetries S.,rn and Cnn are measured. Coo can be determined, up to a fourfold 
ambiguity: 



11 
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1 f / H+- 
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(35) 


- B^%l + 0^)\ ^-[^ 2 ^+ 


±^(g00_^00J(^00^^_^00)(l_O| 
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One obtains the two other solutions Cqq by using 0_ instead of 0+ in the above expres- 
sion. For tto = 90° some simphfications occur: 



0± = ±Jl-Cl^-S^^^ (36) 



^0^ - ^00(1 + o±) { " ^- (^^° " ^"^^ " ^'\ 



±S^^^{B^o _ BZ^){B^X± - ^°°)| . (37) 

This muhivalued function is drawn on Fig. 2, for BABAR, Belle* and the world average, 
with a = 90, 96 and 110°. The B^^ axis spans the range up to 15 x 10~^ despite the 
known experimental upper bound to allow a wide view on the Coo(i3°°) function. Each of 
the two closed loops corresponds to a given choice for the sign of cos(2Q;eflf) cos(2ao): the 
plain lines correspond to 0+, the dotted ones to 0_. The upper and lower parts of each 
closed loop correspond to the ± sign of Eq.(35). The vertical lines indicate the upper and 
lower bounds of the allowed B^^ range. The outermost vertical lines correspond to the 
reBound lower (and upper) limit, while the inner lines correspond to the tough Bound 
limits. 

For Belle*, C^^ + S*^^ being taken equal to unity and then the reBound being equal to 
the tough Bound, the two closed loops are superimposed. 

For BABAR and the world average, both the closed loops and their upper and lower 
components are widely split apart. For BABAR and the W.A., the inner part open when a 
becomes larger, it is the contrary for the outer curves: the larger a becomes, the closer 
they get from each other. 

The variation of the value of a affects the Cqo{B^^) function domain for BABAR, Belle 
and W.A. In the experimentally allowed region {B^^ < 4 x 10"^), both negative and 
positive values are possible for Coo- Here again, such indications are to be taken with 
a grain of salt since the uncertainties on the measurements are not accounted for in the 
discussion (but they are in Section 6). 

4.3 Explicit expression for a 

The identity 2a = 2acs + + leads to the following expression for a: 

sin(2acff)c + cos(2aeflf)s + s ,„„, 

tanQ; = — ^-— — (38) 

cos(2acff)c — sin[2acs)s + c 
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Figure 2: The multivalued function Cqq = /(iS"*"", i3+°, i3°°, S'tttt, C'tttt, «) is drawn as a 
function of B^^ using a = 90 (top), 96 (middle) as measured by the standard CKM fit, 
and 110° (bottom), for the measurements from BABAR (left). Belle* (middle), and the 
world average (right). The outer limits correspond to cos(2a;cff) cos(2ao) > (i.e. the 
reBound of Eq.(24)) and the inner limits to cos(2a;cff) cos(2q;o) < (i.e. the tough Bound 
ofEq.(31)). 



where all quantities on the right hand side can be expressed in term of the observables as 
follow: 



sin(2acff) 
cos(2acff) 



Ott 



1-^L 



±y 1 - sin^(2acflf) 

,,, g+o + g+-(l + a,)/2 - ffl°(l + Coo) 

COS(0) = 

2B+-B+Hl + C^^ 



(39) 

(40) 
(41) 
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sin(0) = ±Vl - c2 (42) 

cos(0) ^ ^" + ^"7^-^-)/^-^"(^-^°°) (43) 



s = sin(0) = ±Vl - c2 (44) 

where the eightfold ambiguity in the range [0, vr] is exphcited by the three arbitrary signs. 
One recovers the GLSS bound by requesting that both c and c are smaller than one. At 
the edge of the allowed B^^ range, for i3°° = iSQLgg.i, one gets: 

«cff = S (45) 



= 


= = 0- 


> a 


Coo = 


1 

»00 
^GLSSil 


(B 


»oo 


" "GLSSiI • 





-a. + ^/S+0S+-(l-Vl-C2^jJ (46) 

(47) 

Being in an Universe where B^^ sits right at the GLSS bound may be viewed, at first 
glance, as good news. Indeed, in that case one has the nice equality a = acs, as if there 
were no Penguins, and, as an odd bonus, one still savors direct CP violation since neither 
Ctttt nor Coo has to vanish. However, to be in an Universe where i3°° is simply close to the 
GLSS bound is bad news! Indeed, in that case, in the neighbourhood of a ^ acs lurk the 
three mirror solutions (2)-(3)-(4), while the other four solutions (5)-(6)-(7)-(8) are thrown 
apart. Furthermore, by construction, when one is close to the GLSS bound dB^^/d6 = 0, 
therefore the locations of the mirror solutions (2) and (4) vary wildly with mild changes 
of i3°°. This subject is discussed in details in Section 5.1 and illustrated on Fig. 3. The 
implication that B^^ = -Bglssi forces a = acs is reminiscent of the Grossman-Quinn 
bound, which is discussed in the next section. 

4.4 The Grossman-Quinn bound revisited 

The Grossman-Quinn (GQ) bound was introduced assuming that no knowledge of a is 
available and Coo is not measured. It is mostly known (and used) as: 

^00 

sm\a - a,s) < ^ • (48) 

This bound appears quite peculiar owing to the results of the previous section where it was 
shown that a equals a^s if '8''° reaches the GLSS bound, not if i3°° is small. In effect, the 
GLSS bound is not necessarily a small number: one may thus be able to set a tight limit on 
the difference | a — acs \ even though B^^ is not small. Mathematically, the constraint on 
I a — acs I does not stem from the size of the branching ratio, but from the proximity of the 
latter to its lower (or upper) bound, if no knowledge is assumed on a. Physically, there is 
no reason to expect a ~ acs except if Penguin diagrams turn out to be very small, which is 
not a strong hope. If it is not the case then only an unfathomable numerical accident may 
lead to i3°° ~ -Sglss;!- ^^ ^^ict, a straight maximization of sin^(2(Q; — aes)) = sin^(0 + 0) 
with respect to Coo yields: 
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C„„|max| = -a, {-^\ g,„^g|^_g„„ I (49) 



/kOO _ »00 \ /kOO _ »00\ 

sm [a - ttcflf j < '^GQ = , (50) 

2B+-B+\l-Cl^ 



This is not a new result: although the end product appears in a quite different form, it 
is derived in the same way in Ref.[19], and is approached in Ref.[18]. It is also present, 
albeit on a low key and in a very different looking form, in the original paper [17] itself ^. 
To recover the GQ bound Eq.(48) from Eq.(50), one may write: 




goo 




{b^^ + ^\- b+^Y 


B+0 


2B+^B+-Jl-Cl^ 



^GQ = ^+ ,. _ -1 UtTT-o - ......., . J • (51) 



With BABAR, Belle and W.A. values both (competing) additionnal terms are small by 
chance. The second term is very small for two reasons: because (for BABMi) C^^ ~ 0.04, 
and because (for all) B^^/B~^^ ~ 0.4 ~ 1/2. The third term is small because i3°° is closed 
to i3+o - B+-/2. 

Hence the validity of the approximation of Eq.(48), when it applies, also deprives the 
GQ bound of its interest: on the one hand, it yields a weak bound, and on the other 
hand, the bound is irrelevant since B^^ and Coo are likely to be measurable, making the 
full isospin analysis tractable. 

Notwithstanding the above considerations, to use Eq.(48) entails an error which can 
go in either directions, as can be seen from Eq.(51). It can lead to overestimate^ or to 
underestimate the bound on | a — acs |. Although with the present data, the effect on 
the calculated values (Table 3) is, by chance, not a large one, there is no reason not to 
use the Sqq bound: the underlying assumptions are identical. 

Before concluding this section, it is appropriate to remind that setting a bound on 
I Oi — cteff I is not an alternative to the full isospin analysis: the latter encompasses the 
first. It is also worth repeating that a tight bound on | a — acs \ also implies that three 
mirror solutions are zeroing on a: it is good news if they all get much closer than the 
statistical power of the isospin analysis, it is very bad news if they just stay at bay within 
one or two sigma's of a. In that case, the measurement is spoilt. 

Remark: Finally, one should keep in mind that Eq.(50) does not refer to ctQ but to 
any of the eight solutions for a. In particular, since ao and S are degenerate in the isospin 
analysis, 6 obeys the constraint of Eq.(50) as well. More precisely, if one considers for 
definiteness the first four mirror solutions for a, and if for example solution (2) reaches 



^Except that, in the second parenthesis, ,8"° is replaced by its lower bound /BgLgg-i, which underesti- 
mates the strength of the bound on | a — a^s \ . 

^The Charles bound [18] and the second GQ bound [17] are intermediate results between the GQ and 
the GLSS bounds which can only understiniate the strenght of the bound on | a — ctcff | • 
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the bound 5'gq, then: 



sm (a(^ij - aeflfj - •^gq- , -^qq ^^^ -^^ ;5oo ^ y^^) 



1 - yi - cL / (ggiss:u - g°°) + (g°° - ggj 



while the locations of the two other mirror solutions are fixed by Eqs.(17) and (18). In 
particular, if C,r7r is small, omitting the last parenthesis which is close to unity, the drift 
from tteflf of a(l) is very small: 

sin2(a(l)-aeff) = 5GQ% ■ (53) 

However, this small value is not a bound. It is only a snapshot: the value reached by a(l) 
when a{2) hits the Sgq bound. If tt(2) is not at the bound sin^(Q;(l) — acflf) takes higher 
or lower values. In particular, a symmetrical situation is encountered when a(l) hits 5'gq. 

5 Perspective for the a measurement 

To illustrate the above relations, different scenarii are considered with possible values of 
B^^, Coo and parameters errors in the future. In Section 5.1 the evolution of the eight 
mirror solutions are studied as a function of i3°°. In Section 5.2, the present world average 
values are used with their errors scaled down by a dreamy factor 10. That is to say a 
luminosity of roughly 10 ab~ is considered both for BABAR and Belle. This deliberately 
unrealistic treatment is chosen to underline the fact that most of the difficulties to predict 
the outcome of the future analyses arise from the vast array of possibilities still open for 
the true characteristics of the B ^ tttt decays. They are not primarily due to the detailed 
performances of the analyses, although of course, the latter are a matter of deep concern. 
In Section 5.3 to illustrate this point, the BABAR values are considered and the present 
analysis performances applied to a 500 fb~ data sample are considered for two theoretical 
scenarii concerning B^^: a nasty one (scenario I) and a nice one (scenario II). 

The results of the analysis are expressed as confidence levels on a. A fit is performed 
minimizing the straight x^ 



X 



, _ (Bt-,-B+-V , (BtZ-B^'\' , (Bl%-B''^ 



001 



to obtain Figs. 4, 5 and 6. 

5.1 Scenario space 

The space of scenarii allowed by the current state of knowledge is a large one. A wide 
point of view is taken and six two-dimensional subsets of the parameters are embraced. 
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These are defined by taking at face value tlie BABAR, Belle* and W.A. measurements, 
together with a = 96° (1.67 rad). 

On Fig. 3 are shown the location of the 8 mirror solutions (Table 2) as a function of 
the value assumed for B^^. Each set of curves starts at the minimal value given by the 
reBound. For the sake of clarity, only the solutions obtained with cos(2Q;eff) negative are 
considered, in agreement with the finding of Section 4.1.2. Otherwise, one should also 
consider another set of six curves, this time starting at the tough Bound. 

Four of the mirror solutions given in Table 2 do not depend on i3°°: beside the true 
solution, they are the solutions (3), (5) and (7). For BABAR as for Belle* one of these 
branching ratio independent solutions is close to the true solution. Whatever B^^, this 
nearby mirror solution has a harmful effect on the measurement. However, the current 
uncertainties on Cyrvr and 3.,^^^ are large enough to allow for this mirror solution to drift 
away swiftly with the coming updates. For the world average, all mirror solutions are 
harmlessly far apart. 

The location of the other four mirror solutions depend on B^^, because Coo is a mul- 
tivalued function of this branching ratio. The Coo{B^^) function drawn on Fig. 2 shows 
that there exist two values of Coo for a given i3°° and for a given sign of the product 
cos(2Q;cfr) cos(2q;o). The three sets of curves on the left hand side are obtained with the 
upper corresponding values of Coo and the three sets of curves on the right hand side with 
the lower values of Coo (for cos(2Q;cff) cos(2q;o) positive). The curves for the two different 
Coo bounce back together when B^'^ hits the reBound B^^^. 

The worst case scenario is when three mirror solutions gather around the true solution 
(ie. when B^° ~ -Bg^Lssd)- With BABAR values, because i3°';\+ ~ -Bg^Lss:! ^r " = 96", 
this happens for i3°° close to the reBound, where one observes that 3 mirror solutions 
are lumped together with the true solution. In that case, the higher i3°°, the better 
the measurement of a. With Belle*, or with the world average values, the worst case 
scenario is avoided because the GLSS bound is well below the reBound. However, nasty 
situations where some mirror solutions collide with a can still occur. With Belle* values, 
this happens by accident (i.e. for -B°° ^ i3g}^ss:i) ^^ -^"^ - 2-4 x 10~^ 



5.2 Guessestimates for high luminosity (10 ab ) 





i3+" (10-6) 


^+0 (10-6) 


^00 (10-6) 


10 ab"' 


4.55 ±0.10 


5.18 ±0.18 


1.9±0.11 


1000 fb"' 


4.55 ±0.19 


5.18 ±0.32 


1.9 ±0.22 




^TTTT 


^TTTT 


Coo 


10 ab"' 


-0.46 ±0.04 


-0.73 ±0.03 


0.8/ -0.38 ±0.08 


1000 fb"' 


-0.46 ±0.07 


-0.73 ±0.07 


0.8/ -0.38 ±0.26 



Table 4: Input values to the fits for the branching ratios, the CP violating asymmetries 
and their extrapolated errors for 10 and 1 af . The central values are the W.A. ones. 
The two values quoted for Coo (i''"^ the ones obtained for cos(2acff) cos(2ao) positive. 
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Figure 3: Location of the 8 mirror solutions as a function of B^^ . The six sets of curves 
refer to BABAR (top), Belle* (middle) and the world average (bottom). Each set starts 
at the minimal value i3°° = B'^y+- Only the solutions obtained with cos(2aoff) negative 
(Section 4-i-3) o^re considered. The curves on the left hand side are obtained with the 
Coo value corresponding to the upper part of the Coo('B°°) function (Fig. 2), the ones on 
the right hand side with the Coo value corresponding to the lower part of the Cqq{B^^) 
function. The vertical lines drawn on each plot represent the true value of a (9Q°) and 
the mirror solutions (3), (5) and (1) which are independent of B'^^ . Four scenarii serve 
as illustrations in Section 5.2 (W.A. using i3°° = 1.9 x 10^^ j and Section 5.3 ('BABAR, 
using B^^ = 1.7 and 0.7 x 10~^^. The corresponding horizontal lines are drawn here to 
indicate the locations of the 8 mirror solutions. 
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Figure 4: Confidence Level CL(a)=Prob{x'^ij^, 1) obtained for W.A. central values at a 
luminosity of 10 a&^ when performing a direct fit with a fine-grain exploration of the 
parameter space (left). The value chosen for Cqo is 0.8. One observes the 8 mirror 
solutions showing up as 8 CL=1 peaks. The right plot is a zoom around the true solution 
a = 96° ('1. 67 rad/ The presence on the left hand side of the nearby mirror solution 
(2) makes the CL curve asymmetrical. The 32% CL is reached at a = 1.67 ± o"± with 
0"+ = 0.07 and a^ = 0.08 (radians). The doted line represents CL(a) for a luminosity of 
1000 fb^^ . In that case, one practically looses all constraints on the left hand side: the 
mirror solution (2) is too close to be resolved. However, the error on the right hand side 
remains useful: cr_ = 0.19. 



The value a = 96° is chosen together with cos(2Q;eflf) negative (Section 4.1.2). The 
central values for W.A. are used as inputs and the errors are extrapolated to 10 ab~^ 
(Table 4). Such a luminosity gives a number of tt^tt" events close to 5000. This set of 
parameters is consistent with two different values of Cqo as shown Fig. 2, which are equal 
to 0.80 and —0.38. The estimated errors for both Cqo are taken to be 0.08. 

The resuh of the fit is shown on Fig. 4 (Cqo = 0.8) and Fig. 5 (Coo = -0.38). Six 
of the mirror solutions appear above a = 96°. The domain allowed at more than 5% 
CL is quite large. Because of that one may be tempted to conclude that tttt does not 
bring much to CKM metrology. However, one should keep in mind that CKM fits receive 
contributions from numerous measurements. It is thus misleading to compare a (tttt) to 
all the other measurements taken at once. Furthermore, precisely because one expects the 
domain in a allowed by a global fit to be quite small in the high luminosity area, the low 
CL values are not anymore relevant. For metrology, what matters then is the sharpness 
of the 8 peaks. If no mirror solution shows up in the vicinity of the correct a value, then 
TTTT may contribute to CKM metrology in a very significant way by helping to pinpoint 
precisely a. This is the case for Cqo = 0.8 where the 8 peaks are sharp, but this is not 
the case for Cqo = —0.38 where the peaks are dull. 

A zoom around a = 96° is shown on Fig. 4 (right) in plain line. Because of the nearby 
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1.6 1.8 

a (rad) 



Figure 5: Confidence Level CL(a)=Prob{Xmini 1) obtained for W.A. central values at a 
luminosity of 10 ab^ when performing a direct fit with a fine-grain exploration of the 
parameter space. The value chosen for Cqq is —0.38. The CL curves obtained in that case 
are to be compared with the ones if Fig. 4 (Cqo = OS). In the present case, the peaks 
merged partly and are much less pronounced. 

mirror solution (2), the CL curve is not symmetrical. Defining a-j- by CL(1.67±o"-|-) = 32% 
one gets a^ = 0.07 and a^ = 0.08. The negative error blows up fast with decreasing 
luminosity, because the fit is quickly not able anymore to resolve the true solution from 
the mirror solution (2). The positive error scales almost with the square root of the 
luminosity because the CL is almost at zero down the pic. The CL(a) for a luminosity 
of 1000 fb~ is shown as well in dotted line. For this case the positive error a^ = 0.19 
is slightly smaller than the (ill-defined) actual precision on a from the global CKM fit 
(= 0.23). Obviously, before then, the measurement of Arn^ by the Tevatron experiments 
will have changed the global CKM picture. 

5.3 Two scenarii at 500 fb^ for BABAR 

Now, only the BABAR central values are considered, at a luminosity of 500 fb~^ which is 
planned to be reached in 2006 by BABAR. Considering the upper left plot on Fig. 3, one 
considers two scenarii for i3°°. The worst case scenarii are obtained when several mirror 
solutions gather around the true solution, staying about one cr[a] away. With BABAR 
values, this happens for B^^ close to the reBound value i3°°^ = 0.6 x 10^^, where 3 mirror 
solutions are lumped together with the true solution. As discussed in Section 4.3 this 
happens because acs being close to a, the reBound is close to the GLSS bound where the 
mirror solutions are exactly degenerate. To show the effect, one considers a scenario with 
the barely acceptable value i3°° = 0.7 x 10^®: this is the scenario I. Scenario II is one 
where, taking the opposite stand, B^^ is far from the reBound. The degeneracy between 
the mirror solutions is then lifted which lets presume that the measurement of a will be 
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B+- (10-6) 


i3+o (10-6) 


^00 (10-6) 


scenario I 


4.7 ±0.26 


5.5 ±0.43 


0.7 ±0.29 


scenario II 


4.7 ±0.26 


5.5 ±0.43 


1.7 ±0.29 




^TTTT 


•^TTTT 


Coo 


scenario I 


-0.19 ±0.09 


-0.40 ±0.10 


0.47 ±1.29 


scenario II 


-0.19 ±0.09 


-0.40 ±0.10 


0.29 ±0.58 



Table 5: Input values to the fits for the branching ratios, the CP violating asymmetries 
and their extrapolated errors. Scenarii I and II use the BABAR values for B^ , B^^, C-n-^r 
and St,t^ and two distinct i3°° guesses. 

easier. The value B^^ = 1.7 x 10-6 is chosen for this case. 

The values used for these two scenarii are presented in Table 5. The result of the fits 
are shown Fig. 6. A zoom of the fits around the true value of a is displayed on the right. 

Whereas, in scenario II, the larger value of the branching ratio B^^ leaves some power 
to the analysis, in scenario I, the vicinity of the mirror solutions deprives the analysis of 
most of its interest. 



6 Standard Model Tests: the (23°°, Cqo) plot 

The above discussion offers a perspective on the future capabilities of tttt analysis to 
constraint the angle a: everything can happen in the few coming years, including a 
precise measurement of a from the B — >■ tttt decays. 

Testing the Standard Model is a different topic. Although this question can be ad- 
dressed by comparing the constraints on a from the standard CKM fit to the ones obtained 
from B ^ 7171 decays, it is a clumsy way to proceed. Indeed, one has to deal with the 
mirror solutions which are irrelevant when one is only willing to know whether or not the 
Standard Model can accomodate the set of six measurements i3"^-, B^^, B^^, C-,,.,,, St,t, 
and Coo knowing a from other measurements. Since the last precise measurements are 
going to be i3°° and Cqo (the latter being still out of reach) it is appropriate to probe the 
Standard Model by considering the (i3°6,Coo) plane. 

The Coo('B*"^) function introduced in Section 4.2 has shown the possible contour, for 
BABAR, Belle* and W.A., of Coo as a function of B^^ (Fig- 2). Here is now presented 
the plane (i3°°,Coo) when a fit is performed using B^~ , B^^, Cyrvr, S-„t, and a, in order 
to evaluate the CL at each point. To look at this plane is a proper way to handle data 
because no bound is needed, nor dealing with mirror solutions. The results for the three 
experimental cases are presented on Fig. 7 for the present luminosity and a = 96" with 
the curves of Fig. 2 superimposed. Considering the allowed region (^3°° < 4 x 10-6)^ one 
observes that there is practically no constraints on Coo, neither for BABAR, Belle* or W.A. 
Similarly, with the present level of knowledge, any value of B^'^ can be accomodated by 
the Standard Model, except exceedingly small ones. 

Figure 8 represent the plane (i3°6. Coo) considering the W.A. results, for a luminosity of 
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1) obtained for scenario I and scenario II. 
The right plot is a zoom around the value 



a 



96° ^1.67 rad; for B^^ = 0.7 x 10"^ (dotted line) and i3"" = 1.7 x 10~^ (plain line). 



500 and 1000 fb ^. To rule out the Standard Model would imply an accurate measurement 
of a large negative value of Cqq and/or a very small value of jB°°. 

7 Conclusion 

The constraints on the angle a are revisited, considering the three branching ratios B^ , 
B^^ and i3°°, and the three CP parameters Cmr, S-^^ and Coo- A new bound is established 
on i3°° taking advantage of the constraints on a available from the standard CKM fit. 
The function Coo('B°°) is introduced and the use of the (i3°°, Coo) plane is advocated as 
a means to probe new physics, when the measurement of Coo will be available. The 
Grossman-Quinn bound is discussed and a compact formula is provided which highlights 
its domain of validity. 

When considering BABAR results, three cases can be considered: 1) if i3°° reaches 
exactly its lower bound, there is no ambiguity for a (but this is academic), 2) if B^^ is 
only close to its lower bound, the locations of the mirror solutions are varying rapidly as 
a function of B^^ and the evaluation of a becomes difficult, 3) if i3°° is far enough from 
the lower bound, the degeneracy between the mirror solutions tend to be lifted and the 
measurement of a may become feasible with only a medium statistics. The current period 
is very exciting, anything can happen in the next three years: including a rather precise 
determination of a. 

The introduction of the (;B°°, Coo) plane in order to test the Standard Model shows 
that it will be difficult to rule out this latter. 
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Figure 7: The (;B°°,Coo) plane is presented for BABAR (left), Belle* (middle) and W.A. 
(right), using the actual luminosity and a = 96°. The curves superimposed are the ones 
given by Eq.(35) for cos(2aeflf) negative (see Fig. 2). 



o 




o 



B''°(10"^) 




Figure 8: The (i3°'^,Coo) plane is presented using the world average branching ratios and 
CP parameters, for 500 fb~^ (left) and 1000 fb^^ (right). 
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